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Abstract. Moments of the characteristic polynomial of a random matrix taken from any of 
the three ensembles, orthogonal, unitary or symplectic, are given either as a determinant 
or a pfaffian or as a sum of determinants. For gaussian ensembles comparing the two 
expressions of the same moment one gets two remarkable identities, one between an n x n 
determinant and an m x m determinant and another between the pfaffian of a 2n x 2n 
anti-symmetric matrix and a sum ofmxm determinants. 



1. Introduction 

Three ensembles of random matrices A have been extensively investigated, namely when 
A is real symmetric, complex hermitian or quaternion self-dual. Recently the behaviour 
of their characteristic polynomial det(xJ — A) has been of interest and we study here its 
integer moments. The same method can be adapted to compute the crossed averages 
(n^det^J-A)). 

Consider a non-negative function w(x) with all its moments finite, 

x m w(x)dx < oo m = 0, 1, ... . (1-1) 

With this weight function w(x) let us define three scalar products, one symmetric and two 
anti-symmetric, as [1] 



(/, 9)2 := J f(x)g(x)w(x)dx = (g, f) 2 (1.2) 
(/, 9)4, ■= J [f(x)g'(x) - f'(x)g{x)]w{x)dx = -{g, /) 4 (1.3) 
= / J f(x)g(y)sign(y - x) v / w(x)w(y)dxdy = -{g,f)i (1.4) 



(f,g) 



and introduce polynomials C n , Q n and R n of degree n, satisfying the orthogonality rela- 
tions 

(C n , C m )2 = c n 5 n ^ m (1-5) 
(Q2m Q2m)4: = (Q2n+1, Q2m+1)4: = (Q2m Q2m+l)4 (1.6) 
(R%niR2m)l = (-R2n+1, R2m+l)l = (i?2n, R2m+l)l = f-n^n,m- \\-*) 
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We will take these polynomials to be monic, i.e. the coefficient of the highest power will 
be taken to be one. The above conditions determine completely the C n (x), Q2n( x ) an d 
R2n{ x )i while Q2n+i( x ) can be replaced by Q2n+i( x ) + a Q2n( x ) with an arbitrary constant 
a. Similarly, R 2n+ i(x) can be replaced by R2n+i(x) + aR 2n (x) with arbitrary a. We will 
choose these constants so that the coefficient of x 2n in Q2n+i( x ) and in R2 n +i(x) is zero, 
thus fixing them also completely. 

The subscript (3 = 1, 2 or 4 is used to remind that it was convenient to use these poly- 
nomials to express the correlation functions of the eigenvalues of real symmetric, complex 
hermitian and quaternion self-dual random matrices respectively [2]. The polynomials 
Cj(x) are said to be orthogonal while the polynomials Qj(x) and Rj(x) are said to be 
skew-orthogonal. 

In what follows n and m are non-negative integers. It can be checked that all the 
given formulas remain valid on replacing a non-existing sum by 0, a non-existing product, 
integral or determinant by 1 and forgetting non-existing rows or columns in a matrix. 

It has long been known that the orthogonal polynomial C n (x) can be expressed as a 
multiple integral [3] 

C n (x)oc I ... j AlwflKx-yjMy^dyj] (1.8) 
J J j=i 

where 

y := {yi,-,y n } (1.9) 

and 

A n (y):= II (Vk-Vi)- ( L1 °) 

l<j<k<n 

Similar multiple integral expressions for the skew-orthgonal polynomials Q n (x) and R n (x) 
were recently discovered by B. Eynard [4] 

Q2n(x) oc I ... j A 4 n (y) f[[(x - yj ) 2 w(yj)dyj} (1.11) 
Q 2n+1 (x)oc j ...j A 4 n (y) lx + 2j2yA f[[(x ~ vtiMv^Vi] i 1 - 12 ) 



2n 



R2n(x)cx J ...j \A 2n (y)\f[[(x - yj )yf^)d yj ] (1.13) 
R2n+l(%) oc J ...J \A 2n (y)\ L + j^yA f[[(x - Vj ) yf^)d yj ]. (1.14) 



2 



It is also known [5] that 



I 2 (n, m; x) :=[... f Aj(y) f[ [(x - yj ) m w( yj )d yj ] (1.15) 
J J j=i 

can be expressed as an m x m determinant. We report here that similar integrals 

Ip(n,m;x) := I ... /" |A n (y)| /3 JJ {( x " 2/j) m [^(^)] 1/(1+5/3,l) ^} (1-16) 

for (3 = 1 or 4, can as well be expressed as a sum ofmxm determinants. They are the 
m-th moments of the characteristic polynomial of an n x n random matrix A 

the parameter (3 taking the values 1, 2 or 4 according as A is real symmetric, complex 
hermitian or quaternion self-dual and the probability density of the eigenvalues y being 
I An(y)P rij l =i[' u; (^j)] 1 ^ 1+ ' 5 ' 3 ' 1 ' ) - The special case w(y) = e~ ay arises when the probability 
density of the random matrix is invariant under any change of basis and the algebraically 
independent parameters specifying the matrix elements are also statistically independent. 

Instead of the real symmetric, complex hermitian and quaternion self-dual matrices 
some authors considered the corresponding ensembles of unitary matrices [6]. The m-th 
moment of the characteristic polynomial in x = exp(ict) 

- / ri \ e%di - ei9k \ p n eiQ - e%Bi \ mde j\ ( Lis ) 

^° l<j<fc<n j=l 

is a constant independent of a which was evaluated using Selberg's integral [6]. In the 
limit of large n with fixed m, this constant has a close relation with the moments of the 
absolute values of the Riemann zeta function ((z) or of the L-functions on the critical line 
Re z — 1/2. Similarly, the moments of the derivative with respect to x of the characteristic 
polynomial are related to the moments of the absolute values of the derivative ('(z) on the 
critical line [6] . 

In case of the real symmetric, complex hermitian and quaternion self-dual matrices, 
which are the only ones we consider here, the m-th moment of the characteristic polynomial 
is a polynomial of order mn, and we do not know whether their zeros are related with the 
Riemann zeta function. 



2. Calculation of I/?(n, m; x) 

The case (3 = 2 is the simplest involving the better known determinants, while the other 
two cases (3 = 1 and 4 are similar involving the less familiar pfaffians. 
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The Ip(n,m;x) can be expressed in two different forms: 

(i) a determinant or a pfaffian of a matrix of size which depends on n; i.e. 

- an n x n determinant in case (3 = 2; 

- a pfaffian of a 2n x 2n anti-symmetric matrix in case (3 = 4; 

- a pfaffian ofanxn (resp. (n + 1) x (n + 1)) anti-symmetric matrix in case /? = 1 for n 
even (resp. odd); 

(ii) a sum of determinants of matrices of size which depends on m; i.e. 

- an m x m determinant in case (3 = 2; 

- a sum ofmxm determinants in case (3 = 4 for m even; 

- a sum ofmxm (resp. (m + 1) x (m + 1)) determinants in case j3 = 1 for n even (resp. 
odd). 

The usefulness of the expressions above, depending on the form of the weight function 
w(x), may be limited. When the weight is gaussian, two of the forms supplying alternate 
expressions for li{n,m\x) as n x n and m x m determinants lead us to an interesting 
known identity. Similarly, the two expressions of I^in, m; x) lead us to an identity between 
a pfaffian of an 2n x 2n anti-symmetric matrix and a sum ofmxm determinants. 



2.1 I/3(n, m; x) as a determinant or a pfaffian of a matrix of size depending on n 

Choose Pj(x) and rcj(x) any monic polynomials of degree j. 
2.1.1 The case 13 = 2 

Write A n (y) as an n x n Vandermonde determinant 

A B (y) = det [y*- 1 ] . fc=1> n = det [P fc -i(%)], )fc=1 ,..., n 

= det[7r fc _ 1 (y i )] i>fc=1 v .. >n . (2.1.1) 



Then 



A* (y) J] [/(Vi)<foi] =/-/ det l p k-i(yj)] det fak-xCy,)] J] [fivMVo 

j = l J J j=l 

= n\ ... / det [irk-iiVj)} II [ P i-i(%)/(%')^] 
J J i=i 

= n! / ... / det [Pj ;_ i ( % - ) 7r fc _ i (yj )/(%•)] dy, 

3=1 



where 



= n\ det p, :hk [ j h „ „ . (2.1.2) 

2;j ,fc := 1 P J (y)n k (y)f(y)dy. (2.1.3) 



In the first three lines of the above equation the indices j, k take the values from 1 to n. 
In the second line we have replaced the first determinant by a single term which is allowed 
by the symmetry of the integrand in the yj; while the last line is obtained by integrating 
over the yj each of them occuring only in one column. 

Setting f(y) = (x — y) m w(y), we get a first form of l2(n, m; x) in terms of an n x n 
determinant 

I 2 (n, m; x) = n\ det [<t>2;j,k{ x )]j, k =o,...,n-i ( 2 - L4 ) 

where 

<h;j,k{x) ■■= J PM*k{y){x - y) m w(y)dy. (2.1.5) 

2.1.2 The case (3 = 4 

Write A^(y) as a 2n x 2n determinant [7] 

At(y)=det[rf, h V ^] j=l x ^ ^ 

= WVi), PL(n)] 1=1 „ ; t =o 2„-i • < 2 - L 6) 

Each variable occurs in two columns, a column of monic polynomials and a column of its 
derivatives. Expand the determinant and integrate to see that [8] 



/ - / A n(y) II fivMVi = ± ^4;ji 1 i 2 04;j3,j4---04;j 2n _ 1 ,i 2 „ (2-1.7) 
J J J = l 

= ««pf[^,*U=o,...,2»-i ( 2 - L8 ) 

where 

<t>4-j,k ■■= J [Pj(y)PL(y) - p;(y)Pk(y)] f(y)dy = -^ k , r (2.1.9) 

In equation (2.1.7) the sum is over all permutations | ? ^ "' ^ n . ^ | of the 2n 

V^l - J2n J 

indices (0, 1, 2n — 1) with j\ < j 2 , 3in-\ < J2n, the sign being plus or minus 
according as this permutation is even or odd. In equation (2.1.8) pf means the pfaffian. 

Setting f(y) = (x — y) m w(y) as before, we get a first form of l4(n,m;x) in terms of 
the pfaffian of a 2n x 2n anti-symmetric matrix 

h(n, m; x) = n\ pf [<t>*;j,k( x )]j,k=o,...,2n-i (2.1.10) 

where 

^,k{x) := j [Pj(y)PL(y) - P;(y)Pk(y)} (* - y) m w(y)d y . (2.1.11) 



2.1.3 The case (3=1 
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The difficulty of the absolute value of A n (y) can be overcome by ordering the variables. 
Writing 



My) II fivi) = det [ p k-i(yj)f(yj)}j, k =i, 



(2.1.12) 



with f(y) = (x — y) m y 'w{y) and using a result of de Bruijn [9], one gets a first form of 
Ii(n,m;x) in terms of the pfaffian of an n x n (resp. (n + 1) x (n + 1)) anti-symmetric 
matrix for n even (resp. odd) 

nl pf [<f>l;3,k( X )]j,k=0,..., n -l 



n even 



ii(n, m; x) 



n\ pf 



<Pl;j,k( X ) a j( X ) 

-a k {x) 



n odd 



(2.1.13) 



j,k=0,...,n-l 



where 



<t>i;j,k{x) := J J P J (z)P k (y)(x - z) m (x - y) m y/w(z)w(y)siga(y - z)dzdy (2.1.14) 



and 



ctj(x) := J P J (y)(x-y) m v ^(yjdy. 



(2.1.15) 



Another way to get the same result is [8] to order the variables, integrate over the 
alternate ones, remove the ordering over the remaining alternate variables and observe 
that the result is given by equations (2.1.13)-(2.1.15). 



2.1.4 An alternative of the expressions above for (3 = 1, 2 and 4 

If one replaces the monic polynomials Pj(y) and nj(y) by Pj(y — x) and iij(y—x), equations 
(2.1.4), (2.1.10) and (2.1.13) become 



h(n,m;x) = (-l) mn n! det [p 2;i) k (x)\ . k=0j n _ 1 
7 4 (n,m; x) = (-l) mn n\ pf [<fi4;j,k( x )]j,k=o,...,2n-i 

nl pf [<Pl;j,k( X )]j,k=0,..., n -l 



Ii(n, m; x) = 



(-l) m n! pf 



<Pi-,j,k( x ) d A x ) 

-e k {x) o 



n even 
n odd 



i,fe=0,...,n-l 



where 



<P2-,j,k(x) := J Pj(y - x)n k (y - x)(y - x) m w{y)dy 



(2.1.16) 
(2.1.17) 

(2.1.18) 



(2.1.19) 



J [P 3 (y - x)P' k (y -x)- P;(y - x)P k (y - x)] (y - x) m w(y)dy (2.1.20) 

<Pi-,j,k( x ) := J J p j( z ~ x ) p k(y ~ x )[{z ~ X )(V - x)] m ^w(z)w(y)sign(y - z)dzdy 

(2.1.21) 

0j(x) := J Pj(y -x)(y- x) m ^y)dy. (2.1.22) 



2.2 I/3 (n, m;x) as determinants of size depending on m 
2.2.1 The case (3 = 2 

For l2(n,m; x), equation (1.15), let us write the integrand as the product of two determi- 
nants: A n (y) given in equation (2.1.1) and 



A n (y)H(x - Vj r = b det 

3=1 



b det 



d 

dx 



x 



k-1 



'Pk-iiy&PHUx) 



j = l,...,n 
k=l,...,n+m 
Z=0,...,m — 1 

j=l,...,n 
k=l,...,n+m 
Z=0,...,m — 1 



(2.2.1) 



where 



f m — 1 



"= n« 



(2.2.2) 



z=o 



Pk(x) is any monic polynomial of degree and pff\x) is its Z-th derivative. Expression 
(2.2.1) is an (m + n) x (m + n) determinant the last m columns of which are Pk{x) and 
its successive derivatives [7]. As the integrand in equation (1.15) is symmetric in the yj, 
we can replace the first n x n determinant by a single term and multiply the result by n\ 



hin, m;x) = b n\ 



det 



,{...{ 



det 



j = l,...,n 
k=l,...,n+m 
1=0,. ..,m — 1 

(0 



II l p 3-i(yj) w (yj) d yj 



Pj-i(yj) p k-i(yj) w (yj), p k-i( x 



j=l,...,n 
k=l,...,n+m 
l=0,...,m — l 



Y[dyj 



(2.2.3) 

As each variable yj occurs only in one column, we can integrate over them independently 



l2(n, m;x)=b n\ det 



J P j _ 1 (y)P k _ 1 (y)w(y)dy 1 Pt l l 1 (x) 



j = l,...,n 
k=l,...,n+m 



(2.2.4) 



If we choose the polynomials Pj(x) to be the orthogonal polynomials Cj(x), equation (1.5), 
then one gets a second form for I 2 {n,m]x) 



J 2 (ra, m;x) =b n\ c ...c n _i det C^ +k (x) 
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k,l=0,...,m — l 



(2.2.5) 



an m x m determinant whose first column consists of C n+ k(x), k = 0, 1, m — 1, and 
the other m — 1 columns are the successive derivatives of the first column. Therefore from 
equation (1.17) 



(det (a;/ - A) m ) = b det 



(2.2.6) 



- k,l=0,...,m — 1 

This result appears in reference [5] as equation (15) apart from a sign which seems to be 
wrong. 



2.2.2 The case /3 = 4 



This method applies only for even moments. For ^(n, 2m; x), equation (1.16), one can 
write the integrand as a single determinant [7] 



A 4 n (y)l[(x-y) 2m = b det 

j=i 



b det 



Vi, ky 



fc-i 



fix 



^(vi),^(vj),n ( V) 



j = l,...,n 
fc=0,...,n+m— 1 
l=0,...,m-l 

j=l,...,n 
k=0,...,n+m—l 
l=0,...,m — l 



(2.2.7) 



where b is given by equation (2.2.2). Each variable yj occurs in two columns. Expanding 
the determinant and integrating one sees that the result is a sum of products of the form 



ifr a s 1 ,s 2 --- a s 2n -i,s 2n det 



>(i-i) 



(x) 



k,l=l,...,m 



(2.2.8) 



with 



a jlk := J [P 3 {y)P' k {y) - P;(y)Pk(y)} w(y)dy = (Pj, P k } 4 (2.2.9) 

the indices s±, S2 n , ji, J m are all distinct, chosen from 0, 1, 2n + m — 1 and the 
sign is plus or minus according as the permutation 



1 

si s 2 



2n-l 2n 

S2n jl 



2n + m — 1 

Jm 



(2.2.10) 



is even or odd. If we choose the polynomials Pj(x) to be the skew-orthogonal polynomials 
Qj(x), equation (1.6), then all the a s j except a2 S ,2s+i = Qs will be zero and we get a 
second form of I^n, 2m; x) in terms of a sum of m x m determinants 



h(n,2m;x) = b n\^q Sl ...q Sn det Qj k X) (x) 



(«) 



k,l=l,...,m 



(2.2.11) 



where the sum is over all choices of s\ < ... < s n such that 2s±, 2s\ + 1, 2s n , 2s n + 1, 
ji, •••) jm are all the indices from to 2n + m — 1 and moreover j\ < ... < j m . 
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2.2.3 The case (3 = 1 



For Ii(n,m;x), equation (1.16), the absolute value sign of A n (y) is the main difficulty. 
Ordering the variables and using the same method as in section 2.2.2 above [7], one has 

h(n,m;x) =/•••/ l A n(y)lH (x - y 3 ) m sjw^y~)dy 3 
J J j=1 L 

= n\ f ... / A n (y) f[ \{x - y 3 ) m ^y~)dy 

J Jy 1 <...<y n j = 1 L 

det[p fc ( % ),P« (x) 



— b n\ 



n 



yi<...<y ri 



w{y 3 )d Vj 



j=l,...,n 
k=0,...,n+m — l 
Z=0,...,m — 1 



(2.2.12) 



with b as in equation (2.2.2). Integrating successively over the alternate variables yi, 2/3, 
1/5, ... and then removing the restriction over the remaining variables y2, 2/4, ... as indicated 
at the end of section 2.1.3 above, one gets for even n 



Ii(n, m;x)=b 



[n/2}\ 



det 



G k (y 2j ),P k (y 2j ),pj!\x) 



i=l,...,[n/2] 
k=0,...,n+m—l 



[n/2] 

n 



w(y 2 j)dy 2j 



with 



/V2 
Pk{yi)V w (yi) d yi- 



(2.2.13) 



(2.2.14) 



In case n is odd, one more column of the numbers gt '■= J Pk{y)\fw{y)dy appears in the 
determinant just after the column of Pk(y n -i)- 

The present situation is exactly as in the case (3 = 4 with each variable occuring in 
two columns. Expanding the determinant one sees that the result contains the expressions 

= (Pj,Pk)i- (2-2.15) 

If we choose Pj(x) = Rj(x), equation (1.7), then writing the result separately for even and 
odd n for clarity, one has 



Ji(2n, m;x) =b (2n)\y~]r 8l ...r 8n det 

(s) 



R 



(i-i) 



jk 



(x) 



k,l=l,...,m 



(2.2.16) 



where the sum is taken over all choices of si < ... < s n such that 2s i, 2s i + 1, 2s n , 
2s n + 1, ji, j m are all the indices from to 2n + m — 1 and moreover ji < ■■■ < jm and 

h(2n+l,m;x) =b (2n+ l)!£V ai ...r s „ det \g jk , (2.2.17) 

^—^ L Jfc Jfc=l,...,m+1; «=l,...,m 

where now the sum is taken over all choices of s± < ... < s n such that 2s±, 2s± + 1, 2s n , 
2s n + l, 

ji-, jrm jm+i are all the indices from to 2n + m and moreover ji < ••• < jm+i- 

Note that equations (1.8), (1.11) and (1.13) are particular cases (m = 1) of equations 
(2.2.5), (2.2.11) and (2.2.16) respectively. If we shift the index of the last row in the right 
hand side of equation (2.2.7) for m = 1 from n to n + 1, i.e. replace the last row 

[P n ( yj ),P^ yj ),P n (x)} 

by the row 

[P n+1 { yj ),P' n+1 { yj ),P n+1 {x)] 

we get the integrand of equation (1.12). Following the procedure which leads to equation 
(2.2.11) we get equation (1.12). Similarly, if we shift the index of the last row from n to 
n + 1 in the right hand side of equation (2.2.12) for m = 1 and follow the procedure leading 
to equation (2.2.16) we get equation (1.14). 



3. Special case of the gaussian weight 

These formulas for Ip{n,m;x) will have little use if one does not know the polynomials 
C n {x), Qn(x) or R n (x). Fortunately one knows them for almost all classical weights. As 

2 

an example we give them here for the gaussian weight w(x) = e~ x over [— oo, oo] in terms 

2 2 

of Hermite polynomials H n (x) := e x (~d/dx) n e~ x . Their verification is straightforward. 
One has 



C n (x) 


= 2- n H n (x) 




(3.1) 




= 2~ n n\ v ^ 




(3.2) 


Q2n(x) 


n | 
3=0 J ' 


Q 2n+1 (x)=2- 2n - 1 H 2n+1 (x) 


(3.3) 


Qn 


= 2- 2n (2n + l)!v^F 




(3.4) 




= 2- 2n H 2n (x) R 2n+ i(x) = 2- 2n [xH 2n (x) - H' 2n {x)\ 


(3.5) 


r n 


= 2 l - 2n {2n)\^. 




(3.6) 



In this case some results can be given other forms. In particular, using the recurrence 
relation 

2xH n (x) = H n+1 (x) + H'Jx) (3.7) 
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or 

xC n {x) = C n+l {x)+ 1 -C' n {x) (3.8) 

one can replace the C^ +k in the determinant (2.2.5) by C n+ k+i(x) so that 

h(n, m;x)=b n\ co...^-^^)^" 1 )/ 2 det [C n+j+k {x\ fc=0>1 ^ m _ t 

2 _„(„-l)/2 rr» -| 
= W 2 ; det[C n+7 - +fc (a;)]., ni (3.9) 

(_ 2 )-m(m-l)/2 YY^jl J lj l k—0,l,...,m—l V V 

Also for any non-negative integer j [10] 

/oo 
(y - x) j e~ v2 dy = v 7 ^ z J Cj(ix) (3.10) 
-oo 

with % = \f—l. Therefore, choosing Pj(y — x) = Wj(y — x) = (y — x) J in equations (2.1.19) 
and (2.1.20), we get 

<P2;j,k(x) = * m+J+k C m+J+k (ix) (3.11) 

<P±,3,k{x) = z m ^' +fc - 1 (k - j)C m+ j + k-l(ix) (3.12) 

and equations (2.1.16) and (2.1.17) then give 

I 2 (n,m;x) = {-i) mn n\ 7r «/2(_i)»(»-i)/2 det [C m+j+k (ix)] . ik=0j ^ (3.13) 
I 4 (n,m;x) = (-i) mn n\ W 2 pf [(A; - ^C^+^x)} . >fc=0> 2n _ x . (3.14) 

Equations (3.9) and (3.13) give a relation almost symmetric in n and m. Writing 
them again 

T (n m- t\ n/2 2 -n(n-l)/2 rr" -| 

h(m,n;ix) { ' W 2 2— (— !)/ 2 n™ -7 ! 

or equivalent ly 

det[C m+i+fc (ix)] iifc=0> ... >n _ 1 (_2)m(m- ^n^i!" 

Equations (3.15) and (3.16) appear in [11] as equations (4.43) and (4.44). 
Equations (2.2.11) and (3.14) give another identity relating the pfaffian 

Pf [(k — j)C2m+j+k-l]j,k=0,...,2n-l 

to a sum ofmxm determinants. 
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Forrester and Witte [11] state that the identity 



(•CO 

J — oo 




2 

(x — yj) m e~ y 3 dyj 


roo 
J — oo 


./_ 0O 0O |A m (y)p/7n i m =1 


2 

(ix — yj) n e~ y i dyj 



const. 



(3-17) 



is valid for any 7. This is not true for 7 ^ 1 as verified below for {n, m} = {n, 1}, {2, 2} 
and {3, 2}. For values of a and 7 such that the integrals below exist, letting 



/ 



d/*(y)/(y) := 



OO /-OO 



00 </ —00 



|A n (y)| 27 n( e " a " 32 ^)/(y) 



(3.18) 

1 v 

(/(y)> := / d/*(y)/(y) - / d/*(y) (3.19) 

^(y) - = a ' n > °; ^(II^ - ^) m ) ( 3 - 2 °) 

A 1 A 1 



3=1 



3 = 1 



where (IljLi^ ~~ 2/j) m ) is a polynomial in x whose coefficients can be evaluated using the 
expressions given in [12]. 

For any n and m = 1, one has 



»/2] 



7 V (2j)! 



n-2j 



(3.21) 



whereas for n = 1 and any m 

X(7, a, 1, m; x) 



J(7, a, 



•X 



[m/2] 

= E 



(2j)l 



a- j x m ~ 2j 



2j J 2 2 ij\ 
= (4a)- m/2 i- m H m (ix^). 

Therefore for any 7, 7', a and n, one verifies that 

J(7,a,n, l;x) = J(7, a, n, 0; x) 



J (7', a/7, 1, ^) J (7', a/7, 1, 0; ix) 
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(3.22) 



(3.23) 



is a constant independent of x. 
For n = m = 2, one has 

j^a^lolx) = { II {X - ^ )2) =x4 + 2 ^ + 2 <^»* 2 + 
= ^ + -(! " 7> 2 + ^2 (1 + 7 + 7 2 ) 

= X(l/7,g/7,2,2;ia:) 24 
J(l/ 7 ,a/ 7 ,2,0;ix)" 1 ' J 

For n = 3 and m = 2, one has 

KM"®-*" 

„6 i q//„.2\ , ,,/„. „. \\ 4 , o//„.2 2\ , ,i/„.2„. „. \\ 2 , /„.2„.2„.2\ 



) 

(3.25) 



x 6 + A(i _ 2 7 )x 4 + A(l - 7 + 3 7 2 )x 2 + ^3 (1 + 3 7 + 5 7 2 ) 



whereas for n = 2 and m = 3 



= x b + 3(2(yf) + 3(y l2/2 ))x 4 + 3(2(yfy 2 ) + 3(y(y$)x z + (y{y 6 2 ) 

(3.26) 



^ 6 + ^(2 - 7)* 4 + £2 ( 3 " 7 + 7 V " £3 (5 + 3 7 + 7 2 ) 



hence 

J(7,g,3, 2;x) _ J(l/7, a/7, 2, 3; ix) 
J(7,a,3,0;x) 1(1/7,0/7, 2,0;^)' 

Selberg's integral also gives the constant [12] 



(3.27) 



J( 7 , a, n, 0; x) = ^^-^(n-^-^n-^-n/i IT Qi±M. ( 3 .28) 

f=i r ( 1 + 7) 



So, what is probably true is 



J( 7 , a, n,m;x) _ ,^ mn J( 7 , a, n, 0; g) 



X(l/7, a/7, m, n; ix) X(l/7, a/7, m, 0; ix) 
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a constant depending on 7, a, n and m. i.e. the gaussians in the numerator and denomi- 
nator of equation (3.17) should have different variances. In particular, 

h(n,m;x) W* 2 -" a n? =0 (2j)! 

Mm.n;**) 1 J 23-/2 pj- r(1 + j/2) • ^.<suj 



4. Conclusion 

The m-th moment of the characteristic polynomial of an n x n random matrix is expressed 
as: an nxn, equations (2.1.4 and 2.1.16), or as an mxm, equation (2.2.5), determinant, for 
the unitary ensemble {j3 = 2); as a pfafhan of a 2n x 2n anti-symmetric matrix, equations 
(2.1.10 and 2.1.17), or a sum ofmxm determinants in case m is even, equation (2.2.11), 
for the symplectic ensemble (/? = 4); as a pfaffian of an n x n (resp. (n + 1) x (n+ 1)) anti- 
symmetric matrix, equations (2.1.13 and 2.1.18), or a sum ofmxm (resp. (m+1) x (m+1)) 
determinants, equation (2.2.16, resp. 2.2.17), for the orthogonal ensemble {j3 = 1) for n 
even (resp. odd). 

2 

In the gaussian case w(y) = e~ y this leads to the remarkable identity almost symmetric 
in m and n, equation (3.15), 

l2(n, m; x) = const. I2(m, n; ix) (4-1) 

or, equation (3.16), 

det [£Wj+fc(aO] j)fc=0 ,...,„_i = const - det [ H n+j+h{ix)] jjk= Q^ m _ 1 (4.2) 

and another identity expressing pf [(k — j)H2 m +j+k(ix)]j k=Q 2n -v equation (3.14), as 
a sum of m x m determinants, equation (2.2.11). Also probably the remarkable general 
identity (3.29) and in particular (3.30) holds. 
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Does this have some deeper significance and/or applications? 
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